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Based on a sample of 227 million BB pairs collected by the BABAR detector at the PEP-II 
asymmetric-energy B Factory at SLAC, we measure the branching fraction B{B'^ tt'^tt") — (1.17± 
0.32 ± 0.10) X 10"®, and the asymmetry C^o^o = -0.12 ± 0.56 ± 0.06. The B° tx'^-k° signal has 
a significance of S.Oct. We also measure B{B+ tt+tt") = (5.8 ± 0.6 ± 0.4) x 10~^ B{B+ 
K+Ti°) = (12.0 ± 0.7 ± 0.6) X 10"^, and the charge asymmetries A^+^a = -0.01 ± 0.10 ± 0.02 and 
-^if +7r0 — 0.06 ± 0.06 ± 0.01. Using isospin relations we find an upper bound on the angle difference 
|a - Qcff i of 35° at the 90% C.L. 

PACS numbers: 13.25.Hw, 12.15.Hh, ll.30.Er 



In the Standard Model (SM) J;he Cabibbo-Kobayashi- 
Maskawa (CKM) matrix Vqqi [jj] describes the charged- 
current coupUngs in the quark sector. The Unitarity Tri- 
angle is a useful representation of relations between CKM 
matrix elements, and measurements of its sides and an- 
gles provide a stringent test of the SM. Following the 
success in measuring the CKM angle /3 an impor- 
tant challenge for the B Factories is the determination of 
the remaining angles. The extraction of the CKM angle 
a = arg [-Vt(jT/t|,/14(jVu*|3] from the time-dependent CP- 
violating asymmetry in the tt+tt" decay mode is 
complicated by the interference of competing amplitudes 
("tree" and "penguin") with different weak phases. The 
difference between a and ckcff , where aes is derived from 
the time-dependent 5° —^ tt+tt^ CP asymmetry, may be 
evaluated using the isospin- related decays i?" vr^vr" 
and _B+ tt+tt" 0. Here and throughout this Let- 
ter, charge conjugate reactions are included implicitly. 
For _B" 7r°7r° the asymmetry may deviate from zero if 
the tree and penguin amplitudes have different weak and 
strong phases. In the SM the decay B^ tt+tt" is gov- 
erned by a pure tree amplitude since penguin diagrams 
cannot contribute to the 1 — 2 final state; as a result no 
charge asymmetry is expected. The B — > Kit system is a 
rich source of information on the understanding of CP vi- 
olation, as has been illustrated by the recent observation 
of direct CP asymmetry in B° K~^Tr~ decays Jsj. Both 
the rate and asymmetry of the _B+ i^+vr" decay may 
be used to extract constraints on penguin contributions 
to the B Kit amplitudes (6|. 

In this Letter, we report a constraint on St,-^ = — a, 
using the measurement of the asymmetry Ct^o^o and 
updated measurements of the branching fractions for 
TT^TT^ and B^ tt+tt" and the charge asym- 
metry At^+t^o. We also measure the branching frac- 
tion for the 5+ K'^tt'^ decay and its charge asym- 
metry Ak+tt°- The asymmetry C^OttO is defined as 
(|Aoo|^ - |Aoo|^)/(|Aoo|^ + |3oo|^), where Aq^ (Aqo) is 
the B'^{B'^) 7r°7r° decay amplitude. For B^ modes, the 
CP-violating charge asymmetry is defined as ^ = (|^P — 
|A|2)/(|A|2 -t- |^|2), where A (A) is the B+ [B') decay 
amplitude. This study is based on 227 x 10^ T{AS) BB 
decays (on- resonance) , collected with the BABAR detec- 
tor. We also use 16fb~^ of data recorded 40MeV below 
the BB production threshold (off-resonance). 

The BABAR detector is described in Ref. 0. The pri- 



mary components used in this analysis are a tracking sys- 
tem consisting of a five-layer silicon vertex tracker (SVT) 
and a 40-layer drift chamber (DCH) surrounded by a 
1.5 T solenoidal magnet, an electromagnetic calorime- 
ter (EMC) comprising 6580 CsI(Tl) crystals, and a ring 
imaging Cherenkov counter (DIRC). 

Candidate 7r° mesons are reconstructed as pairs of 
photons, spatially separated in the EMC, with an in- 
variant mass m-y-y satisfying 110 < m^^ < IGOMeV/c^. 
The mass resolution is 8 MeV/c^ for high energy (above 
2 GeV) 7r°'s 7j. Photon candidates are required to be 
consistent with the expected lateral shower shape, not 
to be matched to a track, and to have a minimum en- 
ergy of 30 MeV. To reduce the background from false vr" 
candidates, the angle 9j between the photon momentum 
vector in the 7r° rest frame and the tt" flight direction 
is required to satisfy |cos6'-y| < 0.95. Candidate tracks 
are required to be within the tracking fiducial volume, 
to originate from the interaction point, to consist of at 
least 12 DCH hits, and to be associated with at least 6 
Cherenkov photons in the DIRC. 

B meson candidates are reconstructed by combining 
a tt" with a charged pion or kaon (/i+) or by combining 
two 7r° mesons. Two variables, used to isolate the 
7r°7r° and B^ /i+7r° signal events, take advantage of 
the kinematic constraints of B mesons produced at the 
T{AS). The first is the beani-energy-substituted mass 
TOES = \/ is/2 + p, • psy/Ef - p|, where {Ei, p,) is the 
four-momentum of the initial e^e~ system, ps is the B 
candidate momentum, both measured in the laboratory 
frame, and -^/i is the e+e" center-of-mass (CM) energy. 
The second variable is AE = EB — y/s/2, where Eb is the 
B candidate energy in the CM frame. The AE resolution 
for signal is approximately 80 MeV for tt^tt", and 

40 MeV for B+ h+n°. 

The primary source of background is e"*"e^ ^ QQ [q ~ 
u, d, s, c) events where a 7r° or from each jet randomly 
combine to mimic a B decay. This jet-like qq background 
is suppressed by requiring that the angle 9s between the 
sphericity axis of the B candidate and that of the remain- 
ing tracks and photons in the event, in the CM frame, 
satisfy I cos 6*3 1 < 0.7(0.8) for B° tt'^tt° {B+ ^ h+n^). 
The other sources of background are B decays to final 
states containing one vector meson and one pseudoscalar 
meson, where one pion is produced almost at rest in the 
B rest frame and the remaining decay products match 
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the kinematics of a B° ^ Tr^Tr" or i3+ /i+7r° decay. 

For the — > 7r°7r° analysis we restrict the mES-Ai? 
plane to the region with m-Es > 5.2GeV/c^ and \AE\ < 
0.4 GeV. For the on-resonance sample we define the sig- 
nal region as the band in the plane with | Ai5| < 0.2 GeV 
and the sideband region as the rest of the plane exclud- 
ing the region which is also populated with B'^ —^ p'^i:^ 
events. The entire plane for the off-resonance data and 
the sideband region for the on-resonance data are kept 
in the fit in order to constrain the qq background param- 
eters. B^ — > ft,+7r° candidates are selected in the region 
with mEs > 5.22GeV/c2 and -0.11 < A£; < 0.15 GeV. 

For B'^ — > 7r°7r° candidates, the other tracks and clus- 
ters in the event are used to determine whether the other 
B meson (Stag) decays as a B^ or B^ (flavor tag). We use 
a multivariate technique to determine the flavor of the 
Btag meson. Events are assigned to one of several mutu- 
ally exclusive categories based on the estimated mistag 
probability and on the source of tagging information. 

The number of signal B candidates is determined with 
an extended, unbinned maximum- likelihood fit. The 
probability density function (PDF) Vi {xj\di) for a sig- 
nal or background hypothesis is the product of PDFs for 
the variables Xj given the set of parameters di . The like- 
lihood function is a product over the N events of the M 
signal and background hypotheses: 

(M \ ^ r 
k=l / j = l Li=l 

For B^ tt^tt" the coefficients Cij are defined as 
Cij = 5(1 — SjAi)ni, where Sj refers to the sign of the 
flavor i;ag of the other B in the event j and is zero for 
untagged events. The fit parameters rii and Ai are the 
number of events and raw asymmetry for i?*^ — > tt'^tt'^ 
signal, B^ p+7r'^ background, and continuum back- 
ground components. The average of branching fraction 
measurements is used to fix n{B^ /O^ti"") to 32 ±6. 
The raw asymmetry for signal is (1 — 2x)(l — 2ll!)Ct^o^o, 
where x = 0.186 ± 0.004 is the neutral B mixing 
probability, and lu is the mistag probability. 

For B~^ /i+7r° the probability coefficients are Cy = 
^(1 — qjAi)ni, where qj is the charge of the track h in the 
event j. The fit parameters rij and Ai are the number of 
events and asymmetry for B^ tt+tt^ and B^ K^ir^ 
signal, continuum, and B background components. The 
B background yields are fixed to the expected number 
of events using the current world averages of branch- 
ing ratios 0|, which are 18 ± 4 for B^ p+7r^ and 
— > p+tt" combined, and 3±1 events for B'^ p^K^ . 
Uncertainties on these numbers are dominated by the un- 
certainty on selection efficiencies, due to the sensitivity 
to the tight requirement in Ai5. 

The variables Xj used for B^ 7r°7r° are ttjesj Ai?, 
and a Fisher discriminant F. The Fisher discrimi- 
nant is an optimized linear combination of ^^'^ 



^jPiCos^6'i, where pi is the momentum and 6i is the 
angle with respect to the thrust axis of the B candi- 
date, both in the CM frame, for all tracks and neutral 
clusters not used to reconstruct the B meson. For both 
the 5° tt'^tt'^ signal and the B~^ /5+7r° background 
the mEs and AE variables are correlated and therefore a 
two-dimensional PDF from a smoothed, simulated distri- 
bution is used. For the continuum background, the mss 
distribution is modeled as a threshold function |0, and 
the AE distribution as a second-order polynomial. The 
PDF for the F variable is modeled as a parametric step 
function (PSF) for all event components. A PSF is a 
variable width binned distribution whose parameters are 
the heights of each bin. The limits of the ten bins F PSF 
are chosen so that each bin contains 10% of the signal 
sample. For B° 7r°7r° and B+ p+Tr° the F PSF pa- 
rameters are correlated with the flavor tagging, and the 
PSF parameters are different for each tagging category. 
Simulated events are used to determine the PSF distri- 
butions for both 5° 7r°7r° and B^ — + p+Tr^. For qq 
background, the F PSF parameters are free in the fit. 

An additional discriminating variable for B~^ —> h^i:^ 
is the Cherenkov angle 6c of the /i+ track. The PDF 
parameters for mESj AE, 6c, and F for the background 
are determined using the data, while the PDFs for signal 
are found from a combination of simulated events and 
data. The tties and AE distributions for qq events are 
treated as in the 7r°7r° case, with parameters al- 

lowed to vary freely in the fit. For the signal, the toes 
and AE distributions are both modeled as a Gaussian 
distribution with a low-side power law tail whose param- 
eters are determined from simulation. The means of the 
Gaussian components are determined from the fit to the 
5+ — > /i+7r° sample and their values used to tune the 7r° 
energy scale in the — > tt^tt^ analysis. The mean of AE 
for the 5+ K^n^ mode is a function of the kaon labo- 
ratory momentum, since a pion mass hypothesis is used. 
The distribution of F is modeled as a Gaussian function 
with an asymmetric variance for the signal, whose pa- 
rameters are obtained from simulation, and as a double 
Gaussian for the continuum background, whose parame- 
ters are determined in the likelihood fit. The difference 
of the measured and expected values of 6c for the pion 
or kaon hypothesis, divided by the uncertainty on 6c, is 
modeled as a double Gaussian function, whose parame- 
ters are obtained from a control sample of kaon and pion 
tracks, from D*+ D^tt+ , D° K-tt+ decays. 

The result of the maximum-likelihood fit for B^ — > 
tt^ttO is n(B° ^ 7r°7rO) = 61 ± 17 (see Tabled, with a 
corresponding statistical significance of 5.2cr. The asym- 
metry is C^o^o = —0.12 ± 0.56. Shown in Fig. ^ are dis- 
tributions of Toes and F, for signal-enriched samples of 
B^ — + tt'^tt^ candidates. The projections contain 25% and 
68% of the signal, 14% and 17% of the p+Tr" background, 
and 2.2% and 4.4% of the continuum background, for 
Toes and F respectively. 
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TABLE I: The results for the modes B° — * tt'V and /i^vr" are summarized. For each mode, the sample size A'', number 

of signal events Ns, total detection efficiency e, branching fraction B, asymmetry A or C^o^o, and the 90% confidence interval 
for the asymmetry are shown. For C^o^o the confidence interval is obtained inferring minimum coverage inside the physical 
region 1]. The first errors are statistical, the second systematic, with the exception of e whose error is purely systematic. 



Mode 


iV 


Ns 


e{%) 


B{10~^) 


Asymmetry 


(90% C.L.) 






8153 


61 ± 17 


23.5 ± 1.4 


1.17 ±0.32 ±0.10 


-0.12 ±0.56 ±0.06 


[-0.88,0.64] 


B+ 




29950 


379 ± 41 


28.7 ± 1.1 


5.8 ±0.6 ±0.4 


-0.01 ±0.10 ±0.02 


[-0.19,0.21] 


B+ ^ 


A'+tt" 


13165 


682 ± 39 


25.0 ± 1.0 


12.0 ±0.7 ±0.6 


0.06 ±0.06 ±0.01 


[-0.06,0.18] 




FIG. 1: Distributions and PDF projections for B" n^Tr". 
Shown are tties (a) and F (b) for candidates that satisfy an 
optimized requirement on the signal probability, based on all 
variables except the one being plotted. PDF projections are 
shown as a dashed line for qq background, a dotted line for B 
background, and a dashed-dotted line for signal. 



FIG. 2: Distributions and PDF projections for B+ -» /i+vr", 
using the method described in the text. For tues (a) the 
signal distributions are combined, while for AE (b) the sig- 
nal B^ tv'^ty'^ (open circles and dashed-dotted curve) and 
B^ K'^tt'^ (solid circles and curve) are shown separately. 
The insets show the combined background components. 



With changes in the analysis technique to measure the 
CP asymmetry, we now find 44 ± 13 signal events in the 
first 123 million BB events, compared to 46 ±13 found in 
Ref. The additional 104 million BB events dataset 
has a signal of 17 ± 11. The signal rates in these two sub- 
sets agree at the 1.3cr level. This result also reflects an 
improved understanding of high energy tt" detection effi- 
ciency. Using a sample of tt*^ mesons from r+ —> Tr+Tr'^i^^ 
decays we apply a 7r° efficiency correction of 0.99 ± 0.03 
to our GEANT simulation, compared to a correction of 
0.88 ± 0.08 applied in Ref. |l3|. 

For _B+ /i+tt" the likelihood fit results are summa- 
rized in Table Using the event-weighting technique 
described in Ref. (lj| we show signal and background 
projections in Fig.|21 For each event, a weight to be sig- 
nal or background is assigned based on a fit performed 
without the specific variable that is plotted. The result- 
ing distributions are normalized to the event yields, and 
are compared to the PDFs used in the fifll fit. 

Systematic uncertainties on the event yields and CP 
asymmetries are evaluated on data control samples, or by 
varying the fixed parameters and refitting the data. In 
order of decreasing importance, the dominant systemat- 
ics on the 5° tt^tt^ branching fraction arise from the 
uncertainty on the AE resolution, the efficiency of the 
tt" reconstruction, and the uncertainty on B background 
event yields. The significance of the — + 7r°7r° signal 
yield, taking systematic effects into account, is 5. Oct. The 
systematic uncertainty on C^o^o is dominated by the un- 
certainties on the B background asymmetry and tagging 



efficiency. 

For B^ /i+7r° the dominant systematic uncertain- 
ties arise from the F signal PDF parameters, selection 
efficiencies, and the AE resolution. Additional system- 
atics arise from uncertainties on the B background event 
yields and particle identification. The systematic uncer- 
tainty on the charge asymmetries is dominated by the 
1% upper limit on the charge bias in the detector [l5(. 

To extract information on i5^^ we use the isospin re- 
lations 3 in conjunction with BABAR measurements of 
C^+^- = -0.09 ± 0.15 ± 0.04 0, the branching frac- 
tion B{B° TT+TT-) = (4.7 ± 0.6 ± 0.2) x 10"^ 0, the 
5° 7r°7r° and _B+ tt+tt" decay rates and the C^o^o 
values reported here. We scan over all values of [(^ttttI and 
calculate a for the decay amplitudes using the method 
described in Ref. The is converted into a con- 

fidence level shown in Fig. |31 from which we derive an 
upper bound on |(5^7r| of 35° at the 90% C.L. 

In summary, we observe 61 ± 17± 5 B'^ tt^tt^ events 
with a significance of 5. Oct including systematic uncer- 
tainties. This corresponds to a branching fraction of 
B{B° -> ttOtt") = (1.17 ± 0.32 ± 0.10) x lO'^, where 
the first error is statistical and the second is systematic. 
We measure the asymmetry C^o^o — — 0.12±0.56±0.06. 
We report branching fractions B{B^ tt+tt'') = (5.8 ± 
0.6 ± 0.4) X 10-s and B{B+ K+n^) = (12.0 ± 0.7 ± 
0.6) X 10"^. The charge asymmetries are At^+t^o — 
-0.01 ± 0.10 ± 0.02 and Ak+^o = 0.06 ± 0.06 ± 0.01; we 
find no evidence for CP violation. In contrast to the re- 
cent measurements of charge asymmetry in B'^ AT+tt" 
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20 40 60 80 
la-a^,,! (degrees) 

FIG. 3: Constraints on the |(5^,r| in terms of confidence level. 
We find an upper bound on \5-„tt\ of 35° at the 90% C.L. 

decays the Ak+ t^o value reported here is compat- 
ible with zero. We use isospin relations on i? ^ tttt 
decay rates and asymmetries to find an upper bound of 
%^\< 35° at the 90% C.L. 
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